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ABSTRACT 

We extend our analyses of the dark matter (DM) distribution in relaxed clusters to the case of Abell 
383, a luminous X-ray cluster at z=0.189 with a dominant central galaxy and numerous strongly-lensed 
features. Following our earlier papers, we combine strong and weak lensing constraints secured with 
Hubble Space Telescope and Subaru imaging with the radial profile of the stellar velocity dispersion 
of the central galaxy, essential for separating the baryonic mass distribution in the cluster core. 
Hydrostatic mass estimates from Chandra X-ray observations further constrain the solution. These 
combined datasets provide nearly continuous constraints extending from 2 kpc to 1.5 Mpc in radius, 
allowing stringent tests of results from recent numerical simulations. Two key improvements in our 
data and its analysis make this the most robust case yet for a shallow slope (3 of the DM density 
profile /?dm oc r - ' 3 on small scales. First, following deep Keck spectroscopy, we have secured the 
stellar velocity dispersion profile to a radius of 26 kpc for the first time in a lensing cluster. Secondly, 
we improve our previous analysis by adopting a triaxial DM distribution and axisymmetric dynamical 
models. We demonstate that in this remarkably well-constrained system, the logarithmic slope of 
the DM density at small radii is ft < 1.0 (95% confidence). An improved treatment of baryonic 
physics is necessary, but possibly insufficient, to reconcile our observations with the recent results of 
high-resolution simulations. 

Subject headings: dark matter — galaxies: clusters: individual (Abell 383) — gravitational lensing: 
strong — gravitational lensing: weak — X-rays: galaxies: clusters — stars: kine- 
matics and dynamics 



1. INTRODUCTION 

The cold dark matter (CDM) paradigm has been re- 
markably successful at predicting the large-scale distri- 
bution of matter in the Universe as well as its observed 
evolution from the earliest epochs to the present day 
(Springcl ct al. 2006). However, much interest has been 
shown in possible discrepancies that remain on small 
scales between its predictions and the available obser- 
vations. A source of continuing puzzlement relates to 
the observed form of the dark matter (DM) profile seen 
in galaxy clusters. 

Numerical simulations predict CDM halos follow a 
self-similar density profile whose three-dimensional (3D) 
form within a scale radius r s approaches a cusp p cx 
with an inner s lope /3 ~ 1 — 1.3 a t the smallest 
resolved scales (e.g.. iNavarro et ail 119961: iGhigna et al.l 
I2000t IDiemand et al.l T20051 ) . Improved resolutio n has 
suggested modest flattening (|Navarro et al.l l2010l ). but 
only at very small radii. This is in contradiction to ob- 
servations based on the combination of strong lensing and 
stellar kinematics which yielded much flatter inner slopes 
of (3 ~ 0.5 for two well-studied clusters (Sa nd et all 
120081) . On the other hand, a steeper = 92 ±0 .04 has 
been inferred in Abell 1703 ([Richard et al.l 120091 ) . possi- 
bly indicating significant scatter in the inner structure of 
clusters. 
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Recently in INewman et al.l (|2009l hereaft er N09) we 
further deve loped the method introduced by ISand et aLI 
(|2004l20"08l ) by incorporating weak lensing constraints on 
the large-scale mass distribution using Subaru imaging 
of Abell 611. This removed a degeneracy between the 
scale radius r s and j3 and led to the first measurement of 
the DM profile across a dynamic range in cluster-centric 
radius (3 kpc to 3.3 Mpc) comparable to that presently 
achieved in simulations. A shallow cusp with /3 < 0.3 
(68%) was derived. 

Here we make two further improvements in our 
methodology and apply these to Abell 383 (z — 0.189). 
This cluster has a regular optical and X-ray morphol- 
ogy, a remarkably low sub-structure fraction, and a 
dominant, near-cir c ular brightest cluster gal axy (BCG) 
( Smit h et all 120051 ISmith fc TavTorl 120081) . ISand et al.1 
(2004) initially stud ied this cluster a ssuming spherical 
symmetry but later (|Sand et al.ll2008l ) undertook a two- 
dimensional (2D) lensing analysis, deriving f3 = 0.45 ± 
0.2. In addition to removing the scale radius degen- 
eracy discussed above, we have considerably extended 
the range of the stellar kinematic constraints via a deep 
Keck spectrum of the BCG, significantly improving our 
knowledge of the mass distribution on < 30 kpc scales. 
Secondly, we use Chandra X-ray data to determine the 
line-of-sight ellipticity in the mass distribution, thereby 
achieving constraints on a 3D model with minimal un- 
certainties arising from projection effects. 

We adopt a cosmology with Da = 0.7, il m = 0.3 and 
Hq = 7O/170 kms -1 Mpc -1 throughout. 

2. OBSERVATIONAL INGREDIENTS 
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Fig. 1. — Strong lensing constraints: HST /WFC2 image in the 
F702W filter with the BCG and other cluster galaxies subtracted 
for clarity. Two families of multiple image systems with known 
spectroscopic redshifts, each comprising two sources, are identified 
in the legend. Critical lines are plotted for both source planes. 

We first discuss in turn the four observational ingredi- 
ents we use to constrain the distribution of dark matter 
and baryons in Abell 383. 

2.1. Strong Lensing 

Figur e [1] shows the multipl y-im aged sources , tabu- 
lated in lSand et all (l200l[2008h and lSmith et~aH (l2005h . 
which comprise the strong lensing constraints: a ra- 
dial/tangential arc at z spC c = 1-01 and a complex system 
in the southwest. In previous analyses the redshift of the 
latter system was unknown, but following Keck/LRIS 
observations in October 2009 we secured a spectroscopic 
redshift z sp c C = 2.55 for images 3C and 4C. The rij=12 
images of n s =4 sources produces 2(r^ — n s ) = 16 con- 
strain ts. Following N09, we use Lenstoo 10 (IJullo et al.l 
in source plane mode for the strong lensing anal- 
ysis. Consistent with earlier work, we adopt an uncer- 
tainty of CTpos = 0'.'5 for the image positions to account 
for systematic modelling uncertainty. 

2.2. Weak Lensing 

The large-scale shear arising from Abell 383 has been 
analyzed with multi-color imaging taken using Suprime- 
Cam at the Subaru telescope. The shear was measured 
in Rc images taken by the authors on 12-13 November 
2007 in excellent seeing of 0'.'57. Broad-band photometry 
from BViz i mages in the S ubaru archive was used with 
the BPZ code (jBenitezI [20001 to obtain photometric red- 
shifts. The procedures closely followed those discussed 
by N09. From a sample of galaxies with 5cr detections in 
Rc, a population with 0.5 < Zf, < 2 was selected for shear 
measurement, yielding a surface density of 25 arcmin -2 . 
As discussed in N09, our shear measurements were cali- 
brated using the recovery fa ctor mwr, = 0.81 ±0 .04 based 
on the STEP2 simulations (jMassev et al.ll2007[ ). 

2.3. Stellar kinematics 

We have subst antially improved the data used by 
ISand et al.l ((2008) by securing a much deeper spectro- 
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scopic exposure of the BCG (Figure H]). Earlier data 
yielded a stellar velocity dispersion profile o~i os {R) ex- 
tending to R = 5 kpc in 3 spatial bins (Figure 4 of 
ISand et alJl2008D . The present data comprise a 22.8 ks 
integration with Keck/LRIS taken on 12 October 2009 
using the 600 mm -1 blue grism and the 600 mm -1 red 
grating blazed at 7500 A. The l'/5 slit yielded a resolution 
of (Tinst = 153 km s _1 at the G band. 

The improved depth of this exposure has enabled us 
to secure a reliable dispersion profile to 26 kpc (cir- 
cularized), which can be verified independently using 
both the G band and Fe absorption lines (Figure [2]). 
Spectra of G and K giants from the MILES library 
(San chez-Blazquez et all 120061) w ere used to synthesiz e 
the optimal stellar template (see iNewman eFaT1[2010l) . 
The gain over earlier data is substantial both in extent 
and sampling. Although ai os (R) h as been measured t o 
very large radii in local cD galaxies iKel -on el al.| [2002: , 
this is to our knowledge the most extended measurement 
yet obtained in a lensing cluster. 

2.4. X-ray 

The final ingredient is the mass distribution probed 
by X-ray emissi on from the i ntrac luster medium. This 
was measured by I Allen et al.l ()2008l hereafter A08) using 
Chandra data. Although their analysis assumed spheri- 
cal symmetry, this has very little effect on the inferred 
spherically-averaged mass profile, as we discuss further 
in Section [5] Mock observations of simulated clusters 
show that non-thermal sources of pressure c ause X-ray- 
deriv e d masses to be biased by ~ —10% (Naga fet al.l 
120071 lLau et al.l 12009ft . To account for this, we place 
a Gaussian prior on mx = MnsE/^truc = 0.9 ± 0.1, 
where Mhse represents the A08 results and M truc the 
true spherically- averaged mass distribution. The 10% 
uncertainty in mx reflects the cluster-to-cluster scatter in 
non-thermal support, as well as uncertainty in th e Chan- 
dra temperature calibration ()Reese eFaT1[20l0l ). From 
the A08 mass profile, we take five points spaced log- 
uniformly from 50 — 600 kpc to match the number of 
independent temperature measurements. (The results 
are not sensitive to the inner limit.) Finally, we add 10% 
in quadrature to the uncertainty in each data point to al- 
low for systematic errors with radial gradients (e.g., non- 
thermal pressure and errors arising from non-sphericity) . 

3. DERIVING THE MASS DISTRIBUTION 

Our model of the cluster mass distribution comprises 
three components: (1) the cluster-scale halo, (2) stellar 
mass in the BCG, and (3) mass in other cluster galax- 
ies which are incorporated as perturbations in the lens- 
ing analysis. The third component is modelled as de- 
scribed in our previous work, including two individually- 
modelled galaxies (PI and P2 in Figured]). Following 
N09, the cluster halo and BCG are described by general- 
ized NFW (gNFW) and dual pseudoisothermal ellipsoid 
(dPIE) profiles, with a key improvement: the newly- 
incorporated X-ray data allows us to consider triaxial 
mass models. 

By combining X-ray and lensing constraints, we can 
directly measure the elongation of the DM distribution 
along the line of sight (l.o.s), thus addressing a key sys- 
tematic uncertainty in deriving the mass density pro- 
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Fig. 2. — Considerably improved kinematic data for the central galaxy, (a) Long slit configuration for the 22.8 ks LRIS exposure; the 8'.'8 
radial extent over which velocity dispersions could be derived is indicated by the vertical red line, (b) Averaged stellar velocity dispersion 
profile as derived independently from the G band and Fe A5270 absorption line regions (blue and red symbols) on either side of the center 
(diamond and plus symbols). The rapid rise in the velocity dispersion at large radii due to the DM halo is clearly evident, (c) Spectra 
for the inner- and outermost spatial bins around the G band. Red curves show the fits to the broadened stellar template, with residuals 
plotted below. 



file. Whereas lensing probes the projected mass con- 
tained in cylinders (2D), the X-ray data is sensitive 
to the spherically-a veraged (3D) enclosed masses (e.g., 
iMorandi et al.ll2~010f ) . The combination thus provides in- 
formation on the l.o.s. geometry. Since the surface den- 
sity of Abell 383 and the BCG isophotes are both nearly 
circular (b/a ~ 0.9), any large departure from sphericity 
must be along the l.o.s. 

In detail, we adopt a triaxial gNFW form for the DM 
halo: 



Pdm(?"£,dm) — 

where 



1 



Po 



9dm (r/r s y(l + r/r s f-0'- 



(1) 



(2) 



es,DM)x 2 + (1 + ££,DM)y 2 + (z/qDM) 2 



Here the z-axis is the l.o.s.; the factor 1/(?dm in Equation 
[1] therefore ensures that the surface density is constant 
as quM varies. The ellipticity of the mass surface density 
£e,dm is related to that of the lensing potential dm 
following iGolse fe Kneibl (|2002t ). Note that q DM > 1 and 
< 1 correspond to prolate and oblate cases, respectively. 

Following N09, the stellar mass of the B CG is modelled 
by a dPIE profile (lEh'asdottir et al.|[2007h fit to HST sur- 
face photometry However, wc consider a more general 
triaxial deprojection, with r e .» defined as in Equation [21 
replacing £e,dm and (/dm by and q«, respectively. 

Before describing our detailed analysis, it is useful 
to gain some physical insight into the effects of vary- 
ing (/dm and q*. qdm governs the ratio between 2D 



and 3D halo masses and is therefore well-constrained 
by the combination of lensing and X-ray data. In con- 
trast, q* is not constrained by long-slit kinematic data, 
and we therefore adopt a prior distribution based on 
kno wledge of the intrinsic axis ratios of elliptical galax- 
ies (jTremblav fe Merrittl 119951 ). We can expect that 
1 < q* < quu (in the prolate notation), both because 
sim ulated DM halos a re much more flattened ((c/a) ~ 
0.5, Uing fe Sutoll2002D than stars in ell iptical galaxies 
((c/a) ~ 0.7. iTremblav fe Merrittl Il995f ). and because 
isotropic dissipation processes in the baryon-dominated 
regime should yi eld rounder mass distributions (e.g., 
lAbadi et al.ll2010l ). Qualitatively, we expect that for a 
fixed halo, rounder stellar orbits will enclose less mass, 
thereby reducing the observed stellar velocity dispersion. 
This introduces a degeneracy between and /3. By ac- 
counting for this degeneracy, we incorporate uncertain- 
ties arising from triaxiality and projection into our final 
results. 

Technically, models are compared to the X-ray data by 
computing spherically-enclosed masses in a triaxial mass 
distribution, as justified in Sections 12 .41 and [51 Fully tri- 
axial dynamical models are not computationally feasible. 
However, since the observations imply DM axis ratios of 
x : y : z ~ 1 : 1.1 : 2, a spheroidal treatment with the 
symmetry axis along the l.o.s. is a very good approxi- 
mation. This represents a significant improvement over 
our previous spherical dynamical models. By assuming a 
two-integral distribution function f (E,L?}, the dy nam- 
ics can be compu ted as described bv lGavazz i (2005|) and 
iQian et all (|1995h . 

Table 1 summarizes the model parameters and our as- 



4 



Newman et al. 



sumed priors. As detailed in N09, models are proposed 
by Markov Chain Monte Carlo (MCMC), and their like- 
lihood is computed as the product of the likelihoods of 
the four datasets. 

4. RESULTS: A TRIAXIAL MODEL WITH SHALLOW 
INNER SLOPE 

Our model fits are plotted in Figure [3] and summa- 
rized in Table 1. To reconcile the observed velocity dis- 
persions with the lensing and X-ray data, shallow DM 
slopes are preferred. As expected, the formal uncertain- 
ties are increased wit h respect to prev ious models that 
neglected triaxiality (jSand et al.l 120081 ). yet we still ob- 
tain B < 0.70 (68% confidence, < 1.0 at 95%). Inter- 
estingly, the DM halo is found to be elongated along 
the l.o. s., with qp M ~ 2, typical of simulated halos 
ping fc Sutol 120021 ) and consistent with the orientation 
selection bias expected for strong lenses. Furthermore, 
the stellar mass-to-light ratio is found to be in agreement 
with the values inferred from stellar population synthesis 
models, assuming a Salpeter IMF (Treu et al. 2010). 

The parameter degeneracies are illustrated in Figure 01 
It is instructive to see how they have been minimized 
by the combination of observational tools unique to our 
method. By including only lensing constraints, for ex- 
ample, we would obtain weaker constraints on 8 due to 
the unconstrained stellar mass. The degeneracy with r s 
is reduced using weak lensing and X-ray probes at large 
radii. Finally, the DM l.o.s. ellipticity (/dm is determined 
by the combination of X-ray and lensing data. 

Equally important to the inferred model parameters is 
the "goodness of fit." As Figure [3] clearly shows, a rela- 
tively simple model, with a single DM halo characterized 
by a simple functional form, fits all the data remarkably 
well, including the extended velocity dispersion data and 
the detailed strong lensing features. (The best-fitting 
models have image plane rms errors of 0'.'3.) The ve- 
locity dispersion profile is particularly encouraging: its 
shape and normalization are well-matched, which was 
typically challeng ing using spherical dynamical models 
(|Sand et al.ll2008L N09). 

5. REMAINING UNCERTAINTIES 

In this Letter we have incorporated additional data 
and more sophisticated models to address the impact of 
projection. Here we consider the residual systematic un- 
certainties. First, we repeated the analysis with all o~i os 
measurements shifted by 10% to account for systematic 
measurement uncertainties. We note that the mild ra- 
dial orbital anisotropy typically observed at lower red- 
shift can only decrease B (see N09). 

Second, we assumed that the BCG is coincident with 
the center of the DM halo, consistent with the < 3 kpc 
projected offset required by the lensing. A similar 3D 
offset has little effect on the enclosed mass outside ~ 
6 kpc, so we can evaluate the effects of this assumption 
by excluding data within 6 kpc. 

A spatially-constant M*/L was assumed, consistent 
with our non-detection o f a color grad i ent o utside the 
central ~ 1". Following iKelson et al.l ((2002) , we esti- 
mate limits on AAf»/i from those on A(B — R). We 
translate this to an uncertainty on R e and repeat the 
analysis shifting R e . 

Finally, we recall that the A08 X-ray measurements 



assumed sphericity, whereas our mass models are non- 
spherical. By calculating the gas emission in a non- 
spherical halo with (/dm = 2, we estimate that spheri- 
cally deprojecting the X-ray observables biases the in- 
ferred (spherically-averaged) mass profile by only ~ 7% 
typically, consiste nt with previous studies ( Gavazzil l2005l : 
iNagai et al.ll2007T) . As discussed in Section l2T4l this small 
bias is comparable to other systematic uncertainties in- 
herent to X-ray analyses and is within our adopted cal- 
ibration uncertainty. We estimate the impact on our re- 
sults by shifting the X-ray masses accordingly. 

In all cases, the limits on 8 shifted by < 0.06. We con- 
clude that the remaining known systematic uncertain- 
ties are much smaller than the projection uncertainty 
addressed in this Letter. 

6. DISCUSSION 

The shallow inner DM slope we find in Abell 383 is 
difficult to reconcile with results from numerical simula- 
tions. DM-only jV-body sim ulations predict 3 ~ 1 — 1.3. 
(e.g., iDiemand et al.l [2005D . Although recent galaxy- 
scale si mulations haye suggested continuous, mo dest flat- 
tening (jStadel et al.l 120091 : iNavarro et all I2010D . this is 
unlikely to affect our comparison, since the slope be- 
comes shallower than NFW only on very small scales 
(< 0.015r s w 4 kpc) that encompass only the innermost 
velocity dispersion bin; this datum does not dominate 
our results. 

Cosmological hydrodynamical simulations incorporat- 
ing baryonic physics (cooling and feedback) with high 
resolution find that baryon condensation in the clus- 
ter core steepens the DM slope (an effect termed 
"adiabatic contraction"), thereby increasing f3 and 
worsening the dis crepancy with observations (e.g., 
IGnedin et al.l l200i 8 ~ 2: I D uffy et al.l f20ToL 8 ~ 
1.5: iSommer-Larsen fc Limousin! I2009L 8 ~ 1.1). The 
amount of steepening depends on the details of the sub- 
grid implementation, but the sign of the effect is consis- 
tent. 

Interestingly, recent hydrodynamical simulations at 
the galaxy and dwarf galaxy scales have shown that 
baryons can soften DM cusps (e.g., iGovernato et al.1 
120091 : iRomano-Diaz et al.ll2009f ); however, the relevant 
processes appear not to scale to galaxy clusters, which 
have deeper potential wells than dwarf galaxies and are 
less baryon-dominated than L* galaxies. Processes that 
have been suggested to counteract adiabatic c ontraction 
in clu ster cores, such as a late dry merging (|Gao et al.l 
I2004D or dynamical friction by infalling bary onic clumps 
( El-Zant et alll20"o"lL l200i INipoti et alJ l2004) , are appar- 
ently subdominant in cosmological hydrodynamical sim- 
ulations and do not lead to a shallow cusp with B < 1. 
Our results on shallow DM cusps in galaxy clusters thus 
seem to require a revision of our understanding of either 
the DM backbone or the most relevant baryonic physics 
for shaping the cluster core. 

7. SUMMARY 

We have extended our previous analyses of Abell 383 
by incorporating weak lensing and X-ray data. Based on 
deep Keck spectroscopic observations, we have also mea- 
sured - for the first time in a lensing cluster - an extended 
velocity dispersion profile that demonstrates a clear rise 
in the outer regions in response to the cluster potential. 
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Fig. 3. — Top: Mass enclosed in cylinders (2D) for the dark and stellar components, as well as the total mass distribution, with tangential 
reduced shear (g+) data inset at the same radial scale. Bottom: Mass enclosed in spheres (3D), with velocity dispersion data inset and 
X-ray constraints overlaid. All bands show 68% confidence regions. The models acceptably fit all constraints from ~ 2 kpc to ~ 1.5 Mpc. 
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Fig. 4. — Joint posterior probability density of /3 with four other model parameters: M„/Ly, r s , 5dMi an d 5*. The potential degeneracies 
are reduced by the inclusion of complementary constraints. Contours show 68% and 95% confidence regions, and red lines indicate the 
mean fj to highlight the degeneracy slopes. Histograms show marginalized pos teriors. Upper limits o n (3 are indicated in the right panel, 
along with the range (shaded) reported for cluster-scale TV-bo dy simulations by Dicmand ct al. (2005). Dotted lines in the left panel show 
the M- t /Ly inferred from Subaru photometry using kcorrect (Blanton & Rowcis 2007), assuming Chabrier and Salpeter IMFs. 
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TABLE 1 

Models Inferred from Strong and Weak Lensing, Kinematic, and X-ray Data 



Quantity 


Units 


Prior 


Marginalized Posterior 


gNFW DM halo 

4,DM 

Position angle (P.A.) 
r s 

"0,DM 



9DM 


deg 
kpc 
km s" 1 


[0,0.15] 
[-10,30] 
log-uniform* 
log-uniform* 
uniform* 
[1.5,2.7] 


O^™ 

10.5+*;? 

112+ 61 
1 -30 

1629+^° 

59+ ' 30 

1 O7+0.28 
-"-■^ -0.16 


BCG stellar mass 
r C ore,fcut (dPIE) 
b/a 
P.A. 

M./iy 
9* 


kpc 

deg 
(M/L v )q 


0.82, 40.5 
0.88 
15 

[1.8] 

1/?* ~ 0.54 ±0.15 


(fixed) 
(fixed) 
(fixed) 

n QK+0.90 

1 30+ 015 
1 - JU -0.13 


Calibration parameters 
m WL 
m x 




0.81 ±0.04 
0.90 ±0.10 


n 7 o+0.03 
U -' S -0.02 

1 01+"' 04 
1 - ul -o.os 


Cluster galaxy perturbers 
°~0,* 
r C ut,* 
r C ut,Pi 

?"cut,P2 


km s — 1 
kpc 
kpc 
kpc 


159 ± 40 
[20, 80] 
[20, 120] 
[20, 80] 


122+ J? 
t 
t 
t 



Note. — Posteriors are summarized using the mode and the 68% confidence interval. These priors were found to be non- 
restrictive following initial tests with broader priors, [a, b] denotes a uniform prior, while /U ± a denotes a Gaussian, which is 
truncated at 1.5<r for calibration parameters. *In practice, independent, uniform priors are placed on linear combinations of 
(log r s , log op, dm, 0) for effic i ency; these are equivalent to broad, uniform priors on logr s , loguo.DM, and fi. The prior on q, is 
a fit to lTremblav fc Merrittl |l995). Formulae to convert oo to other quantities are given in N09. (To,* and r cut ,* denote dPIE 
parameters for an L* cluster galaxy (see N09). Perturber radii marked with f have nearly flat posteriors, which are omitted. 



Dark Matter Distribution in Abell 383 



7 



As a result of these improved datasets, we have refined 
our modelling to consider triaxial mass distributions. We 
demonstrate that all four observational ingredients are 
essential to obtain a complete three-dimensional view 
of Abell 383 over a very wide dynamic range of three 
decades in radius. 

We find that the DM distribution in Abell 383 is clearly 
triaxial, consistent with ACDM numerical simulations. 
However, its DM profile has a shallow density cusp with 
f3 < 1 (95%), a result that appears inconsistent with 
numerical simulations of clusters in a ACDM context at 
moderate significance. Although Abell 383 represents 
only a single well-studied cluster, comparable data are 
currently in hand to extend this analysis to sample of 9 



clusters whose ensemble properties will be presented in 
a forthcoming paper. 
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